Genome-wide maps of highly-similar intrachromosomal repeats that mediate ectopic recombination in three human genome assemblies.

Repeated sequences spread throughout the genome play important roles in shaping the structure of chromosomes and facilitating the generation of new genomic variation. Through a variety of mechanisms, repeats are involved in generating structural rearrangements such as deletions, duplications, inversions, and translocations, which can have the potential to impact human health. Despite their significance, repetitive regions including tandem repeats, transposable elements, segmental duplications, and low-copy repeats remain a challenge to characterize due to technological limitations inherent to many sequencing methodologies. We performed genome-wide analyses and comparisons of direct and inverted repeated sequences in the latest available human genome reference assemblies including GRCh37 and GRCh38 and the most recent telomere-to-telomere alternate assembly (T2T-CHM13). Overall, the composition and distribution of direct and inverted repeats identified remains similar among the three assemblies but we observed an increase in the number of repeated sequences detected in the T2T-CHM13 assembly versus the reference assemblies. As expected, there is an enrichment of repetitive regions in the short arms of acrocentric chromosomes, which had been previously unresolved in the human genome reference assemblies. We cross-referenced the identified repeats with protein-coding genes across the genome to identify those at risk for being involved in genomic disorders. We observed that certain gene categories, such as olfactory receptors and immune response genes, are enriched among those impacted by repeated sequences likely contributing to human diversity and adaptation. Through this analysis, we have produced a catalogue of direct and inversely oriented repeated sequences across the currently three most widely used human genome assemblies. Bioinformatic analyses of these repeats and their contribution to genome architecture can reveal regions that are most susceptible to genomic instability. Understanding how the architectural genomic features of repeat pairs such as their homology, size and distance can lead to complex genomic rearrangement formation can provide further insights into the molecular mechanisms leading to genomic disorders and genome evolution. Author summary: This study focused on the characterization of intrachromosomal repeated sequences in the human genome that can play important roles in shaping chromosome structure and generating new genomic variation in three human genome assemblies. We observed an increase in the number of repeated sequence pairs detected in the most recent telomere-to-telomere alternate assembly (T2T-CHM13) compared to the reference assemblies (GRCh37 and GRCh38). We observed an enrichment of repeats in the T2T-CHM13 acrocentric chromosomes, which had been previously unresolved. Importantly, our study provides a catalogue of direct and inverted repeated sequences across three commonly used human genome assemblies, which can aid in the understanding of genomic architecture instability, evolution, and disorders. Our analyses provide insights into repetitive regions in the human genome that may contribute to complex genomic rearrangements.

Break-induced replication underlies formation of inverted triplications and generates unexpected diversity in haplotype structures.

Background: The duplication-triplication/inverted-duplication (DUP-TRP/INV-DUP) structure is a type of complex genomic rearrangement (CGR) hypothesized to result from replicative repair of DNA due to replication fork collapse. It is often mediated by a pair of inverted low-copy repeats (LCR) followed by iterative template switches resulting in at least two breakpoint junctions in cis . Although it has been identified as an important mutation signature of pathogenicity for genomic disorders and cancer genomes, its architecture remains unresolved and is predicted to display at least four structural variation (SV) haplotypes. Results: Here we studied the genomic architecture of DUP-TRP/INV-DUP by investigating the genomic DNA of 24 patients with neurodevelopmental disorders identified by array comparative genomic hybridization (aCGH) on whom we found evidence for the existence of 4 out of 4 predicted SV haplotypes. Using a combination of short-read genome sequencing (GS), long- read GS, optical genome mapping and StrandSeq the haplotype structure was resolved in 18 samples. This approach refined the point of template switching between inverted LCRs in 4 samples revealing a DNA segment of ∼2.2-5.5 kb of 100% nucleotide similarity. A prediction model was developed to infer the LCR used to mediate the non-allelic homology repair. Conclusions: These data provide experimental evidence supporting the hypothesis that inverted LCRs act as a recombinant substrate in replication-based repair mechanisms. Such inverted repeats are particularly relevant for formation of copy-number associated inversions, including the DUP-TRP/INV-DUP structures. Moreover, this type of CGR can result in multiple conformers which contributes to generate diverse SV haplotypes in susceptible loci .

Genetic convergence of developmental and epileptic encephalopathies and intellectual disability.

AIM: To determine whether genes that cause developmental and epileptic encephalopathies (DEEs) are more commonly implicated in intellectual disability with epilepsy as a comorbid feature than in intellectual disability only. METHOD: We performed targeted resequencing of 18 genes commonly implicated in DEEs in a cohort of 830 patients with intellectual disability (59% male) and 393 patients with DEEs (52% male). RESULTS: We observed a significant enrichment of pathogenic/likely pathogenic variants in patients with epilepsy and intellectual disability (16 out of 159 in seven genes) compared with intellectual disability only (2 out of 671) (p<1.86×10-10 , odds ratio 37.22, 95% confidence interval 8.60-337.0). INTERPRETATION: We identified seven genes that are more likely to cause epilepsy and intellectual disability than intellectual disability only. Conversely, two genes, GRIN2B and SCN2A, can be implicated in intellectual disability without epilepsy; in these instances intellectual disability is not a secondary consequence of ongoing seizures but rather a primary cause. What this paper adds A subset of genes are more commonly implicated in epilepsy than other neurodevelopmental disorders. GRIN2B and SCN2A are implicated in intellectual disability and epilepsy independently.

Pathogenic Variants in CEP85L Cause Sporadic and Familial Posterior Predominant Lissencephaly.

Lissencephaly (LIS), denoting a "smooth brain," is characterized by the absence of normal cerebral convolutions with abnormalities of cortical thickness. Pathogenic variants in over 20 genes are associated with LIS. The majority of posterior predominant LIS is caused by pathogenic variants in LIS1 (also known as PAFAH1B1), although a significant fraction remains without a known genetic etiology. We now implicate CEP85L as an important cause of posterior predominant LIS, identifying 13 individuals with rare, heterozygous CEP85L variants, including 2 families with autosomal dominant inheritance. We show that CEP85L is a centrosome protein localizing to the pericentriolar material, and knockdown of Cep85l causes a neuronal migration defect in mice. LIS1 also localizes to the centrosome, suggesting that this organelle is key to the mechanism of posterior predominant LIS.

MN1 C-terminal truncation syndrome is a novel neurodevelopmental and craniofacial disorder with partial rhombencephalosynapsis.

MN1 encodes a transcriptional co-regulator without homology to other proteins, previously implicated in acute myeloid leukaemia and development of the palate. Large deletions encompassing MN1 have been reported in individuals with variable neurodevelopmental anomalies and non-specific facial features. We identified a cluster of de novo truncating mutations in MN1 in a cohort of 23 individuals with strikingly similar dysmorphic facial features, especially midface hypoplasia, and intellectual disability with severe expressive language delay. Imaging revealed an atypical form of rhombencephalosynapsis, a distinctive brain malformation characterized by partial or complete loss of the cerebellar vermis with fusion of the cerebellar hemispheres, in 8/10 individuals. Rhombencephalosynapsis has no previously known definitive genetic or environmental causes. Other frequent features included perisylvian polymicrogyria, abnormal posterior clinoid processes and persistent trigeminal artery. MN1 is encoded by only two exons. All mutations, including the recurrent variant p.Arg1295* observed in 8/21 probands, fall in the terminal exon or the extreme 3' region of exon 1, and are therefore predicted to result in escape from nonsense-mediated mRNA decay. This was confirmed in fibroblasts from three individuals. We propose that the condition described here, MN1 C-terminal truncation (MCTT) syndrome, is not due to MN1 haploinsufficiency but rather is the result of dominantly acting C-terminally truncated MN1 protein. Our data show that MN1 plays a critical role in human craniofacial and brain development, and opens the door to understanding the biological mechanisms underlying rhombencephalosynapsis.

BRAT1 encephalopathy: a recessive cause of epilepsy of infancy with migrating focal seizures.

Epilepsy of infancy with migrating focal seizures (EIMFS), one of the most severe developmental and epileptic encephalopathy syndromes, is characterized by seizures that migrate from one hemisphere to the other. EIMFS is genetically heterogeneous with 33 genes. We report five patients with EIMFS caused by recessive BRAT1 variants, identified via next generation sequencing. Recessive pathogenic variants in BRAT1 cause the rigidity and multifocal seizure syndrome, lethal neonatal with hypertonia, microcephaly, and intractable multifocal seizures. The epileptology of BRAT1 encephalopathy has not been well described. All five patients were profoundly impaired with seizure onset in the first week of life and focal seizure migration between hemispheres. We show that BRAT1 is an important recessive cause of EIMFS with onset in the first week of life, profound impairment, and early death. Early recognition of this genetic aetiology will inform management and reproductive counselling.

GGC Repeat Expansion and Exon 1 Methylation of XYLT1 Is a Common Pathogenic Variant in Baratela-Scott Syndrome.

Baratela-Scott syndrome (BSS) is a rare, autosomal-recessive disorder characterized by short stature, facial dysmorphisms, developmental delay, and skeletal dysplasia caused by pathogenic variants in XYLT1. We report clinical and molecular investigation of 10 families (12 individuals) with BSS. Standard sequencing methods identified biallelic pathogenic variants in XYLT1 in only two families. Of the remaining cohort, two probands had no variants and six probands had only a single variant, including four with a heterozygous 3.1 Mb 16p13 deletion encompassing XYLT1 and two with a heterozygous truncating variant. Bisulfite sequencing revealed aberrant hypermethylation in exon 1 of XYLT1, always in trans with the sequence variant or deletion when present; both alleles were methylated in those with no identified variant. Expression of the methylated XYLT1 allele was severely reduced in fibroblasts from two probands. Southern blot studies combined with repeat expansion analysis of genome sequence data showed that the hypermethylation is associated with expansion of a GGC repeat in the XYLT1 promoter region that is not present in the reference genome, confirming that BSS is a trinucleotide repeat expansion disorder. The hypermethylated allele accounts for 50% of disease alleles in our cohort and is not present in 130 control subjects. Our study highlights the importance of investigating non-sequence-based alterations, including epigenetic changes, to identify the missing heritability in genetic disorders.

Rhombencephalosynapsis: Fused cerebellum, confused geneticists.

Rhombencephalosynapsis (RES) is a unique cerebellar malformation characterized by fusion of the cerebellar hemispheres with partial or complete absence of a recognizable cerebellar vermis. Subsets of patients also have other brain malformations such as midbrain fusion with aqueductal stenosis, characteristic craniofacial features (prominent forehead, flat midface, hypertelorism, ear abnormalities), and somatic malformations (heart, kidney, spine, and limb defects). Similar to known genetic brain malformations, the RES cerebellar malformation is highly stereotyped, yet no genetic causes have been identified. Here, we outline our current understanding of the genetic basis for RES, discuss limitations, and outline future approaches to identifying the causes of this fascinating brain malformation.

A mutational comparison of adult and adolescent and young adult (AYA) colon cancer.

BACKGROUND: It is possible that the relative lack of progress in treatment outcomes among adolescent and young adult (AYA) patients with cancer is caused by a difference in disease biology compared with the corresponding diseases in younger and older individuals. There is evidence that colon cancer is more aggressive and has a poorer prognosis in AYA patients than in older adult patients. METHODS: To further understand the molecular basis for this difference, whole-exome sequencing was conducted on a cohort of 30 adult, 30 AYA, and 2 pediatric colon cancers. RESULTS: A statistically significant difference in mutational frequency was observed between AYA and adult samples in 43 genes, including ROBO1, MYC binding protein 2 (MYCBP2), breast cancer 2 (early onset) (BRCA2), MAP3K3, MCPH1, RASGRP3, PTCH1, RAD9B, CTNND1, ATM, NF1; KIT, PTEN, and FBXW7. Many of these mutations were nonsynonymous, missense, stop-gain, or frameshift mutations that were damaging. Next, RNA sequencing was performed on a subset of the samples to confirm the mutations identified by exome sequencing. This confirmation study verified the presence of a significantly greater frequency of damaging mutations in AYA compared with adult colon cancers for 5 of the 43 genes (MYCBP2, BRCA2, PHLPP1, TOPORS, and ATR). CONCLUSIONS: The current results provide the rationale for a more comprehensive study with a larger sample set and experimental validation of the functional impact of the identified variants along with their contribution to the biologic and clinical characteristics of AYA colon cancer. Cancer 2018;124:1070-82. © 2017 American Cancer Society.

Plasmid-Based Materials as Multiplex Quality Controls and Calibrators for Clinical Next-Generation Sequencing Assays.

Although next-generation sequencing technologies have been widely adapted for clinical diagnostic applications, an urgent need exists for multianalyte calibrator materials and controls to evaluate the performance of these assays. Control materials will also play a major role in the assessment, development, and selection of appropriate alignment and variant calling pipelines. We report an approach to provide effective multianalyte controls for next-generation sequencing assays, referred to as the control plasmid spiked-in genome (CPSG). Control plasmids that contain approximately 1000 bases of human genomic sequence with a specific mutation of interest positioned near the middle of the insert and a nearby 6-bp molecular barcode were synthesized, linearized, quantitated, and spiked into genomic DNA derived from formalin-fixed, paraffin-embedded-prepared hapmap cell lines at defined copy number ratios. Serial titration experiments demonstrated the CPSGs performed with similar efficiency of variant detection as formalin-fixed, paraffin-embedded cell line genomic DNA. Repetitive analyses of one lot of CPSGs 90 times during 18 months revealed that the reagents were stable with consistent detection of each of the plasmids at similar variant allele frequencies. CPSGs are designed to work across most next-generation sequencing methods, platforms, and data analysis pipelines. CPSGs are robust controls and can be used to evaluate the performance of different next-generation sequencing diagnostic assays, assess data analysis pipelines, and ensure robust assay performance metrics.

Germline De Novo Mutations in GNB1 Cause Severe Neurodevelopmental Disability, Hypotonia, and Seizures.

Whole-exome sequencing of 13 individuals with developmental delay commonly accompanied by abnormal muscle tone and seizures identified de novo missense mutations enriched within a sub-region of GNB1, a gene encoding the guanine nucleotide-binding protein subunit beta-1, Gß. These 13 individuals were identified among a base of 5,855 individuals recruited for various undiagnosed genetic disorders. The probability of observing 13 or more de novo mutations by chance among 5,855 individuals is very low (p = 7.1 × 10(-21)), implicating GNB1 as a genome-wide-significant disease-associated gene. The majority of these 13 mutations affect known Gß binding sites, which suggests that a likely disease mechanism is through the disruption of the protein interface required for Gα-Gßγ interaction (resulting in a constitutively active Gßγ) or through the disruption of residues relevant for interaction between Gßγ and certain downstream effectors (resulting in reduced interaction with the effectors). Strikingly, 8 of the 13 individuals recruited here for a neurodevelopmental disorder have a germline de novo GNB1 mutation that overlaps a set of five recurrent somatic tumor mutations for which recent functional studies demonstrated a gain-of-function effect due to constitutive activation of G protein downstream signaling cascades for some of the affected residues.

Analytical Validation and Application of a Targeted Next-Generation Sequencing Mutation-Detection Assay for Use in Treatment Assignment in the NCI-MPACT Trial.

Robust and analytically validated assays are essential for clinical studies. We outline an analytical validation study of a targeted next-generation sequencing mutation-detection assay used for patient selection in the National Cancer Institute Molecular Profiling-Based Assignment of Cancer Therapy (NCI-MPACT) trial (NCT01827384). Using DNA samples from normal or tumor cell lines and xenografts with known variants, we assessed the sensitivity, specificity, and reproducibility of the NCI-MPACT assay in five variant types: single-nucleotide variants (SNVs), SNVs at homopolymeric (HP) regions (≥3 identical bases), small insertions/deletions (indels), large indels (gap ≥4 bp), and indels at HP regions. The assay achieved sensitivities of 100% for 64 SNVs, nine SNVs at HP regions, and 11 large indels, 83.33% for six indels, and 93.33% for 15 indels at HP regions. Zero false positives (100% specificity) were found in 380 actionable mutation loci in 96 runs of haplotype map cells. Reproducibility analysis showed 96.3% to 100% intraoperator and 98.1% to 100% interoperator mean concordance in detected variants and 100% reproducibility in treatment selection. To date, 38 tumors have been screened, 34 passed preanalytical quality control, and 18 had actionable mutations for treatment assignment. The NCI-MPACT assay is well suited for its intended investigational use and can serve as a template for developing next-generation sequencing assays for other cancer clinical trial applications.

Using mixtures of biological samples as process controls for RNA-sequencing experiments.

BACKGROUND: Genome-scale "-omics" measurements are challenging to benchmark due to the enormous variety of unique biological molecules involved. Mixtures of previously-characterized samples can be used to benchmark repeatability and reproducibility using component proportions as truth for the measurement. We describe and evaluate experiments characterizing the performance of RNA-sequencing (RNA-Seq) measurements, and discuss cases where mixtures can serve as effective process controls. RESULTS: We apply a linear model to total RNA mixture samples in RNA-seq experiments. This model provides a context for performance benchmarking. The parameters of the model fit to experimental results can be evaluated to assess bias and variability of the measurement of a mixture. A linear model describes the behavior of mixture expression measures and provides a context for performance benchmarking. Residuals from fitting the model to experimental data can be used as a metric for evaluating the effect that an individual step in an experimental process has on the linear response function and precision of the underlying measurement while identifying signals affected by interference from other sources. Effective benchmarking requires well-defined mixtures, which for RNA-Seq requires knowledge of the post-enrichment 'target RNA' content of the individual total RNA components. We demonstrate and evaluate an experimental method suitable for use in genome-scale process control and lay out a method utilizing spike-in controls to determine enriched RNA content of total RNA in samples. CONCLUSIONS: Genome-scale process controls can be derived from mixtures. These controls relate prior knowledge of individual components to a complex mixture, allowing assessment of measurement performance. The target RNA fraction accounts for differential selection of RNA out of variable total RNA samples. Spike-in controls can be utilized to measure this relationship between target RNA content and input total RNA. Our mixture analysis method also enables estimation of the proportions of an unknown mixture, even when component-specific markers are not previously known, whenever pure components are measured alongside the mixture.

Robustness of Next Generation Sequencing on Older Formalin-Fixed Paraffin-Embedded Tissue.

Next Generation Sequencing (NGS) technologies are used to detect somatic mutations in tumors and study germ line variation. Most NGS studies use DNA isolated from whole blood or fresh frozen tissue. However, formalin-fixed paraffin-embedded (FFPE) tissues are one of the most widely available clinical specimens. Their potential utility as a source of DNA for NGS would greatly enhance population-based cancer studies. While preliminary studies suggest FFPE tissue may be used for NGS, the feasibility of using archived FFPE specimens in population based studies and the effect of storage time on these specimens needs to be determined. We conducted a study to determine whether DNA in archived FFPE high-grade ovarian serous adenocarcinomas from Surveillance, Epidemiology and End Results (SEER) registries Residual Tissue Repositories (RTR) was present in sufficient quantity and quality for NGS assays. Fifty-nine FFPE tissues, stored from 3 to 32 years, were obtained from three SEER RTR sites. DNA was extracted, quantified, quality assessed, and subjected to whole exome sequencing (WES). Following DNA extraction, 58 of 59 specimens (98%) yielded DNA and moved on to the library generation step followed by WES. Specimens stored for longer periods of time had significantly lower coverage of the target region (6% lower per 10 years, 95% CI: 3-10%) and lower average read depth (40x lower per 10 years, 95% CI: 18-60), although sufficient quality and quantity of WES data was obtained for data mining. Overall, 90% (53/59) of specimens provided usable NGS data regardless of storage time. This feasibility study demonstrates FFPE specimens acquired from SEER registries after varying lengths of storage time and under varying storage conditions are a promising source of DNA for NGS.

Randomized Trial of Oral Cyclophosphamide and Veliparib in High-Grade Serous Ovarian, Primary Peritoneal, or Fallopian Tube Cancers, or BRCA-Mutant Ovarian Cancer.

PURPOSE: Veliparib, a PARP inhibitor, demonstrated clinical activity in combination with oral cyclophosphamide in patients with BRCA-mutant solid tumors in a phase I trial. To define the relative contribution of PARP inhibition to the observed clinical activity, we conducted a randomized phase II trial to determine the response rate of veliparib in combination with cyclophosphamide compared with cyclophosphamide alone in patients with pretreated BRCA-mutant ovarian cancer or in patients with pretreated primary peritoneal, fallopian tube, or high-grade serous ovarian cancers (HGSOC). EXPERIMENTAL DESIGN: Adult patients were randomized to receive cyclophosphamide alone (50 mg orally once daily) or with veliparib (60 mg orally once daily) in 21-day cycles. Crossover to the combination was allowed at disease progression. RESULTS: Seventy-five patients were enrolled and 72 were evaluable for response; 38 received cyclophosphamide alone and 37 the combination as their initial treatment regimen. Treatment was well tolerated. One complete response was observed in each arm, with three partial responses (PR) in the combination arm and six PRs in the cyclophosphamide alone arm. Genetic sequence and expression analyses were performed for 211 genes involved in DNA repair; none of the detected genetic alterations were significantly associated with treatment benefit. CONCLUSION: This is the first trial that evaluated single-agent, low-dose cyclophosphamide in HGSOC, peritoneal, fallopian tube, and BRCA-mutant ovarian cancers. It was well tolerated and clinical activity was observed; the addition of veliparib at 60 mg daily did not improve either the response rate or the median progression-free survival.

Regulation and expression of the ATP-binding cassette transporter ABCG2 in human embryonic stem cells.

The expression and function of several multidrug transporters (including ABCB1 and ABCG2) have been studied in human cancer cells and in mouse and human adult stem cells. However, the expression of ABCG2 in human embryonic stem cells (hESCs) remains unclear. Limited and contradictory results in the literature from two research groups have raised questions regarding its expression and function. In this study, we used quantitative real-time PCR, Northern blots, whole genome RNA sequencing, Western blots, and immunofluorescence microscopy to study ABCG2 expression in hESCs. We found that full-length ABCG2 mRNA transcripts are expressed in undifferentiated hESC lines. However, ABCG2 protein was undetectable even under embryoid body differentiation or cytotoxic drug induction. Moreover, surface ABCG2 protein was coexpressed with the differentiation marker stage-specific embryonic antigen-1 of hESCs, following constant BMP-4 signaling at days 4 and 6. This expression was tightly correlated with the downregulation of two microRNAs (miRNAs) (i.e., hsa-miR-519c and hsa-miR-520h). Transfection of miRNA mimics and inhibitors of these two miRNAs confirmed their direct involvement in the regulation ABCG2 translation. Our findings clarify the controversy regarding the expression of the ABCG2 gene and also provide new insights into translational control of the expression of membrane transporter mRNAs by miRNAs in hESCs.

Reduced seed region-based off-target activity with lentivirus-mediated RNAi.

Along with silencing intended target genes, transfected siRNAs regulate numerous unintended transcripts through a mechanism in which the equivalent of a microRNA-like seed region in the siRNA recognizes complementary sequences in transcript 3' UTRs. Amelioration of this off-target silencing would lead to more accurate interpretation of RNA interference (RNAi) experiments and thus greatly enhance their value. We tested whether lentivirus-mediated delivery of shRNA is prone to the sequence-based off-target activity prevalent in siRNA experiments. We compared target gene silencing and overall impact on global gene expression caused by multiple sequences delivered as both transfected siRNAs and lentivirus vector-expressed shRNAs. At equivalent levels of target gene silencing, signatures induced by shRNAs were significantly smaller than those induced by cognate siRNAs and arose less frequently from seed region activity. Interestingly, the low level of seed region-based off-target activity exhibited by shRNAs resulted in down-regulation of transcripts that were largely distinct from those regulated by siRNAs. On the basis of these observations, we recommend lentivirus-mediated RNAi for pathway profiling experiments that measure whole genome transcriptional readouts as well as for large-scale screens when resources for extensive follow up are limited.